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Confexf. Oscillatory reconnection is a time-dependent magnetic reconnection mechanism that naturally produces periodic outputs 
from aperiodic drivers. 

Aims. This paper aims to quantify and measure the periodic nature of oscillatory reconnection for the first time. 
Methods. We solve the compressible, resistive, nonlinear MHD equations using 2.5D numerical simulations. 

Results. We identify two distinct periodic regimes: the impulsive and stationary phases. In the impulsive phase, we find the greater 
the amplitude of the initial velocity driver, the longer the resultant current sheet and the earlier its formation. In the stationary phase, 
we find that the oscillations are exponentially decaying and for driving amplitudes 6.3 - 126.2 km/s, we measure stationary-phase 
periods in the range 56.3 - 78.9 s, i.e. these are high frequency (0.01 - 0.02 Hz) oscillations. In both phases, we find that the greater 
the amplitude of the initial velocity driver, the shorter the resultant period, but note that different physical processes and periods are 
associated with both phases. 

Conclusions. We conclude that the oscillatory reconnection mechanism behaves akin to a damped harmonic oscillator. 
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1. Introduction 



Magnetohydrodynamics (MHD) - Waves - Sun: corona - Sun: magnetic topology - Sun: os- 



Traditionally, magnetic reconnection and MHD wave theory 
have been viewed as separate areas of solar physics (see, 
e.g., Priest & Forbes l20U0l Rob erts l2004l De Moortel 120031 
Nakariakov & Verwichte I2005I De Moortel & Nakariakov 
I2012l l. However, this is a misconception: we know that (steady- 
state) reconnection models not only generate outflows/waves, 
but also require inflows/waves (e.g. Parker [19571 Sweet [1958; 
Petschek 1964). Several authors have already challenged this 
point-of-view (e.g. Craig & McClymont 1 199 1 1 Longcope & 
Priest HHl Murray et al. 120091 McLaughlin et al. 120091 12TJT21 
and their investigations contribute to our understanding of dy- 
namic or time-dependent models of magnetic reconnection. Of 
particular importance to this paper is the work of McLaughlin et 
al. ([2009) which is the first demonstration of reconnection nat- 
urally driven by MHD wave propagation, via a process entitled 
oscillatory reconnection. 

MHD wave propagation in inhomogeneous media is a funda- 
mental plasma process and the study of MHD waves in the 
neighbourhood of magnetic null points directly contributes to 
this area (see review by McLaughlin et al. 1201 II ). It is known 
that null points - weaknesses in the magnetic field where the 
field strength, and hence the Alfven speed, is zero - and sepa- 
ratrices - topological features that separate regions of different 
magnetic flux connectivity - are an inevitable consequence of 
the distributed isolated magnetic flux sources at the photospheric 
surface, where the number of such null points will depend upon 
the magnetic complexity of the photospheric flux distribution 
(see, e.g., review by Longcope [2005 , and Close et al. 2004; 
Regnier et al. 2008] Lon gcope & Parnell 2009 for the statistics 
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of coronal null points). It is also now known that MHD wave 
perturbations are omnipresent in the corona (e.g. Tomczyk et al. 
2007). Thus, these two areas of scientific study; MHD waves 
and magnetic topology, will encounter each other in the corona, 
i.e. MHD waves will propagate into the neighbourhood of coro- 
nal null points (e.g. blast waves from a flare will at some point 
encounter a null point). 

McLaughlin & Hood d2004l [20031 12006al I2006bl > investigated 
the behaviour of linear MHD waves (fast and slow magnetoa- 
coustic waves and Alfven waves) in the neighbourhood of a 
variety of 2D null points. It was found that the (linear) fast 
wave is focused towards the null point by a refraction effect 
and all the wave energy, and thus current density, accumulates 
close to the null, i.e. null points will be locations for prefer- 
ential heating by (linear) fast waves. The Alfven wave prop- 
agates along magnetic fieldlines and so accumulates along the 
separatrices (in 2D) or along the spine or fan-plane (in 3D). 
Waves in the neighbourhood of a single 2D null point have also 
been investigated using cylindrical models, in which the gener- 
ated waves encircled the null point (e.g. Bulanov & Syrovatskii 
H9801 C raig & McClvmont [T99Tl [T9931 Craig & Watson [19921 
Hassam 1992) and it was found that the wave propagation leads 
to an exponentially-large increase in the current density (see 
also Ofman [19921 Ofman et al. [19931 Steinolfso n et a l. [19951 
and a comprehensive review by McLaughlin et al. 1201 ll for fur- 
ther details). 3D MHD wave activity about coronal null points 
has been investigated by various authors (e.g. Galsgaard et al. 
120031 Pontin & Galsgaard[2007J P ontin et al .120071 McLaughlin 
et al. 120081 Galsg aard & Pontin 1201 lal 1201 lbl Thurgood & 
McLaughlin [20T21 

Reconnection can occur when strong currents cause the mag- 
netic fieldlines to diffuse through the plasma and change 
their connectivity (Parker 119571 Sweet 119581 Petschek 11964) . 
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However, these papers did not include the effect of gas pres- 
sure, which would act to limit the growth of the current den- 
sity. In considering the relaxation of a 2D X-type neutral point 
disturbed from equilibrium, Craig & McClymont (1 1 99 li t found 
that free magnetic energy is dissipated by a physical mechanism 
which couples resistive diffusion at the null to global advection 
of the outer field, which they called oscillatory reconnection. 
An example of oscillatory reconnection generated by magnetic 
flux emerging into a coronal hole was reported by Murray et 
al. ([2009) who found a series of "reconnection reversals" take 
places as the system searches for equilibrium, i.e. the system 
demonstrates oscillatory reconnection in a self-consistent man- 
ner. The physics behind oscillatory reconnection has been inves- 
tigated by McLaughlin et al. ( 120091 ), Murray et al. ( 120091 ) and 
Threlfall et al. d20T2b . 

McLaughlin et al. ( 120121) investigated the long-term evolu- 
tion of an initially-buoyant magnetic flux tube emerging into a 
gravitationally-stratified coronal hole environment and reported 
on the resulting oscillations and outflows. They found that the 
physical mechanism of oscillatory reconnection naturally gen- 
erates quasi-periodic vertical outflows with a transverse/swaying 
aspect. There is currently a great deal of interest in observations 
of transverse motions in the solar atmosphere (e.g. Tomczyk 
et al. 120071 De Po ntieu e t al. 120071 120TT1 Cirtain et al . [20071 
Er delyi & Taroya n [20081 Nishizuka et al. [20081 I20TT1 He et 
al. l2UU9al I2009bt Liu et al. 120091 I20TT1 Mc intosh et al. I20TT1 
Oka moto & De Pontieu I20TT1 Morton et al. 120121 Yurchyshyn 
et al. 120121 ) and transverse/swaying motions have been observed 
over a range of wavelengths, speeds, temperatures and scales. 
However, the origin of these propagating, transverse oscillations 
remains a mystery, and these authors often note that the chal- 
lenge remains to understand how and where these waves are 
generated in the solar atmosphere. Liu et al. ( 1201 II ) summarises 
possible generation mechanisms for these transverse motions, 
including an oscillating wake from a coronal mass ejection or 
periodic reconnection (see, e.g., Chen & Priest 2006; Sych et 
al. 12009). The physical mechanism of oscillatory reconnection 
is another possible source of these transverse motions. As re- 
ported by McLaughlin et al. (120121 ). the transverse behaviour 
seen in the periodic jets originating from the reconnection re- 
gion of the inverted Y-shaped structure is specifically due to the 
oscillatory reconnection mechanism, and would be absent for a 
single, steady-state reconnection jet. The physical mechanism 
also naturally generates periodic outputs even though no peri- 
odic driver is imposed on the system. 

Thus, there is a clear interest in furthering our understanding of 
the periodic nature of oscillatory reconnection. In this paper, we 
investigate the periodic signal generated by the mechanism, with 
a specific interest in measuring periods and decay rates as well 
as the robustness of our results, i.e. how do the results vary with 
the strength of the driver. 

1.1. Overview of McLaughlin et al. (2009) 

This paper will closely follow the work of McLaughlin et al. 
(2009 ) as we investigate the periodic nature of oscillatory recon- 
nection. These authors investigated the behaviour of nonlinear 
fast magnetoacoustic waves near a 2D X-type neutral point and 
found that the incoming wave deforms the null point into a cusp- 
like point which in turn collapses to a current sheet. The system 
then evolves periodically through a series of horizontal/vertical 
current sheets with associated changes in connectivity, i.e. the 
system demonstrates the mechanism of oscillatory reconnection. 



More specifically, McLaughlin et al. (2009) found that the in- 
coming (fast) wave propagates across the magnetic fieldlines and 
the initial profile, an annulus, contracts as the wave approaches 
the null point. This is the refraction behaviour that is typical 
of fast wave behaviour around magnetic null points (see, e.g., 
McLaughlin et al. I201 II ) and results from the spatially-varying 
(equilibrium) Alfven-speed profile. 

The incoming wave was observed to develop discontinuities (for 
a physical explanation, see Appendix B of McLaughlin et al. 
2009 or, alternatively, Gruszecki et al. l201 1) and these disconti- 
nuities form fast oblique magnetic shock waves, where the shock 
makes B refract away from the normal. Interestingly, the shock 
locally heats the initially /3 = plasma, creating /3 + at these 
locations. 

At a later time, the shocks overlap, forming a shock-cusp, which 
leads to the development of hot jets and in turn these jets sub- 
stantially heat the local plasma and significantly deform the lo- 
cal magnetic field. By the time the shock waves reach the null, 
the (initially X-point) magnetic field has been deformed such 
that the separatrices now touch one another rather than intersect- 
ing at a non-zero angle (Priest & Cowley [19751 call this 'cusp- 
like'). The osculating field structure continues to collapse and 
forms a horizontal current sheet. However, the separatrices con- 
tinue to evolve: the jets at the ends of the (horizontal) current 
sheet continue to heat the local plasma, which in turn expands. 
This expansion squashes and shortens the current sheet, forc- 
ing the separatrices apart. The (squashed) current sheet thus re- 
turns to a 'cusp-like' null point that, due to the continuing ex- 
pansion of the heated plasma, in turn forms a vertical current 
sheet. In effect, the (net) restoring force acts to return the (de- 
formed) null point to its equilibrium state, but overshoots the 
equilibrium. The phenomenon then repeats itself: jets heat the 
plasma at the ends of this newly-formed (vertical) current sheet, 
the local plasma expands, the (vertical) current sheet is short- 
ened, the system attempts to return itself to equilibrium, over- 
shoots and forms a (second) horizontal current sheet. The evolu- 
tion proceeds through a series of horizontal and vertical current 
sheets, and the system clearly displays oscillatory behaviour. It is 
also interesting to note that the final state is non-potential, where 
this is because the plasma to the left and right of the null is (lo- 
cally) hotter than that above and below. Consequently, a thermal- 
pressure gradient exists and causes the X-point to be slightly 
closed up in the vertical direction (i.e. generating a small, pos- 
itive current). It is important to note that the non-potential final 
state is still in force balance and will eventually return to a po- 
tential state, but on a far greater timescale than our simulations 
(^diffusion ~ Rm = 10 4 Alfven times, where R m is the magnetic 
Reynolds number). We also note that there is nothing unique 
about the orientation of the first current sheet being horizontal 
followed by a vertical, this simply results from the particular 
choice of initial condition, and McLaughlin et al. ( 120121 ) use the 
more general terminology orientation 1 and orientation 2. 
McLaughlin et al. (2009 ) also present evidence of reconnection 
occurring in the system; reporting both a change in fieldline con- 
nectivity (qualitative evidence) and changes in the vector po- 
tential which directly showed a cyclic increase and decrease in 
magnetic flux on either side of the separatrices (see their Figures 
12 and 13). Hence, since the system displayed both oscillatory 
behaviour and reconnection, it was concluded that the system 
demonstrated the phenomenon of oscillatory reconnection. 
Our paper has the following outline: the basic setup, equations 
and assumptions are described in ^2) the periodic nature of os- 
cillatory reconnection is detailed in |3] and the conclusions are 
given in ^4] 
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2. Basic Equations 

We consider the nonlinear, compressible, resistive MHD equa- 
tions: 



5v / m " 
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m 

dt 

| + (T-V), 
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= V x (v x B) + 77V 2 B , 
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= - P v ■ v + - ijr + gshock , 

cr 



dp 
dt 



+ V-(pv) = 



(1) 



where p is the mass density, v is the plasma velocity, B the 
magnetic induction (usually called the magnetic field), p is the 
plasma pressure, p = An x 10~ 7 Hm _I is the magnetic permeabil- 
ity, cr is the electrical conductivity, 77 = 1 1 per is the magnetic 
diffusivity, e = p/p(y- 1) is the specific internal energy den- 
sity, y — 5/3 is the ratio of specific heats and j = V x B//x is the 
electric current density. 

We solve these governing equations numerically using a 
Lagrangian remap, shock-capturing code called LARE2D (Arber 
et al. 120011) . which utilizes artificial shock viscosity to introduce 
dissipation at steep gradients. The details of this technique, often 
called Wilkins viscosity, can be found in Wilkins ( 119801 1. Thus, 
6shock represents the viscous heating at shocks. 
We now introduce a change of scale to non-dimensionalise all 
variables. Letting v = V(>v*, B = BW, x = Lx* , y = Ly*, z = Lz* , 
p = pop*, p = pop", j = 70 j*, V = V*/L, t = t f and 77 = 770, 
where * denotes a dimensionless quantity and Vo, B, L, po, po, jo, 
to and 770 are constants with the dimensions of the variable they 
are scaling. We then set Bj yjppo = vq and vo = L/to (this sets vo 
as a constant background Alfven speed). We also set jo = B/pL 
and rjoto/L 2 = /J" 1 , where R m is the magnetic Reynolds number, 
and choose R,„ = 10 4 . This process non-dimensionalises equa- 
tions (Q} and under these scalings, t* = 1 (for example) refers to 
t = to — L/vo; i.e. the time taken to travel a distance L at the 
background Alfven speed. 

There is no fixed dimensional length scale to our X-point sys- 
tem (X-points are scale-free), and thus we have a great deal 
of freedom in choosing our dimensional constants. We choose 
L = 1 Mm and B = 1 G (for simplicity, and where these choices 
allow an intuitive understanding of equation [2] below) and we 
chose a coronal density of po = 5 x 10 _13 kg/m 3 and coronal tem- 
perature of Tq = 10 6 K. This sets vo = BJ -\fppo = 126.2 km/s, 
t = L/vo = 7.93 seconds and j = B/pL = 8 x 10~ 5 A. The 
values returned from equations ([TJ are made dimensional using 
these solar constants. 



2. 1 . Basic equilibrium and numerical set-up 

To set-up our system, we follow the numerical framework of 
McLaughlin et al. (2009). Thus, we consider a simple 2D X- 
type neutral point as our equilibrium magnetic field, where the 
initial field is taken as: 



B = -(y,x,0) 



(2) 



where B — 1 G is a characteristic field strength and L — 1 Mm 
is the length scale for magnetic field variations. This magnetic 
field can be seen in McLaughlin et al. (2009, their Figure 1). 
Initially, we consider the equilibrium plasma to be cold: T — K 
(i.e. B(t = 0) = 0) and, hence, ignore plasma pressure effects. 



However, McLaughlin et al. (2009) showed that magnetic shocks 
will heat the plasma and so the plasma will not remain cold (see 
e.g. §1.5 in Priest & Forbes l20U0l l. 

In order to excite a pure fast magnetoacoustic wave, we consider 
an initial condition that perturbs velocity purely across the equi- 
librium magnetic field. Using the terminology of McLaughlin 
et al. (120091 1. there are three distinguishing velocity components 
that can be considered in our system: 

• v_l = v x Bo ■ z = v x B y - v y B x , which corresponds to the 
velocity across the equilibrium field, and hence corresponds 
to the fast magnetoacoustic wave (the only MHD wave that 
can cross fieldlines), 

• V|| = v • Bo = v x B x + v v By, which corresponds to the velocity 
along the equilibrium field and corresponds to the propaga- 
tion of the slow magnetoacoustic wave, 

• v z = v • t, which corresponds to the velocity in the invariant 
direction and hence corresponds to the Alfven wave. 

Interestingly, these three velocity components - each isolating 
an individual MHD mode - are in good agreement with those 
reported by Thurgood & McLaughlin (2012) who used the equi- 
librium magnetic field and the flux function (which is parallel 
to the invariant direction) to define an orthogonal coordinate 
system to isolate and identify the propagation of each of the 
MHD modes. Using their convention, v_i_ = v • Bo x A, where 
A = j (y 2 - x 2 } z is the flux function, and where perturbations in 

the Bo x A-direction were shown to correspond to those of the 
fast wave (for helicity-free systems, such as our 2D X-point). 
To perturb our system, we consider an initial condition in veloc- 
ity such that: 

v ± (x,y,t = 0) = 2Csin[ff(r-4.5)] for 4.5 < r < 5.5 , 

v,| (x,y, t = 0) = v, (x,y, t = 0) = (3) 

where r 2 — x 2 + y 2 and 2C is our initial amplitude and initial 
condition ([3]) describes a circular, sinusoidal pulse in v ± . Thus, 
as argued above, we initially generate (only) a pure fast wave in 
our system. Note that the velocity profile prescribed by equation 
(01 appears as the symmetric m = mode in v ± , but corresponds 
to the asymmetric m = 2 mode in Cartesian components. This is 
why the first current sheet has horizontal orientation, as per 91.11 
This initial pulse will naturally split into two waves - an outgoing 
wave and an incoming wave - each of amplitude C. In this paper, 
we will focus on the incoming wave, i.e. the wave propagating 
towards the null point. In this paper, we will conduct a param- 
eter study of initial wave amplitude C. Note that setting C — 1 
recovers the results of McLaughlin et al. (2009) and choosing 
a small value for C, say C = 0.001, recovers the the linear re- 
sults from McLaughlin & Hood (120041; i.e. see Appendix A of 
McLaughlin et al. 2009). Under our dimensionalisation, a choice 
of C = 1 corresponds to an incoming wave with maximum ini- 
tial Cartesian velocity v ± /r = vo/5 = 25.2 km/s at a distance of 
r = 5L = 5 Mm and an equilibrium magnetic field strength of 
5G. 

The governing equations (Q3 with initial conditions ([3]) are solved 
computationally in a square domain x, y e [-20,20] with a nu- 
merical resolution of 6144 x 6144. Zero gradient boundary con- 
ditions are applied to the variables B, p, e at the four boundaries, 
and v is set to zero on all boundaries, i.e. reflective boundaries. 
A numerical damping region exists for x 2 + y 2 > 6 which gradu- 
ally removes kinetic energy from the outgoing waves and so all 
oscillations that enter this region are slowly damped away, and 
hence they do not influence the behaviour about the null. The 
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Fig. 1. Plot of time evolution of j z (0, 0, f) (measured in milliAmps) for < t < 480 seconds. Red line indicates impulsive phase and 
black line indicates stationary phase. Insert shows the same time evolution over 158.9 < t < 480 seconds (same horizontal axis, 
different vertical axis). The black dotted lines indicate the formation times of all the horizontal current sheets and the green line 
indicates j'fi na i = 0.8 165 jo = 6.5 x 10~ 5 A. The blue dashed lines indicate an exponentially-damped envelope max(J z )\ stationary x e~ M + 
j'flnai and minC/JIstationary x e~ M + j'flna] where A = -0.015 s" 1 . 



(equilibrium) Alfven speed increases with distance from the null 
point and, hence, waves accelerate as they propagate outwards. 
Since we do not want reflected waves to influence our null point, 
implementation of such a damping region is essential. 

3. Aperiodic driver leading to periodic behaviour 

We set C = 1 in equations 0]) and, as expected, we recover the 
results of McLaughlin et al. (2009) and readers are directed to 
that paper for full details (primarily their Figures 2, 5 & 6) and 
also this paper's 91.11 One of the key results from McLaughlin 
et al. (12009b was the production of a periodic response re- 
sulting from an aperiodic input, i.e. the physical mechanism 
of oscillatory reconnection naturally gave rise to periodic be- 
haviour. This periodic response can be quantitatively measured 
by analysing the time evolution of the (electric) current density, 
specifically j = (0, 0, j z ), at the null point itself. This can be seen 
in McLaughlin et al. (2009, their Figure 10) and the analysis 
of this time signal is the primary focus of this paper. In recre- 
ating the simulations of McLaughlin et al. (2009J, we recover 
this periodic time series, but at a higher numerical resolution, 
6144 x 6144, and at a higher cadence, dt - 0.05 f ( >. The time 
evolution of j z (0, 0, t) for < t < 480 s can be seen in Figure[T] 
Note that due to the symmetry of our system, the null is always 
located at the origin. 

We identify two distinct regimes in Figure Q] < t < 90 s (red 
line in Figure Q]) which we refer to as the impulsive phase and 
t > 90 s (black line in Figure[TJ which we refer to as the station- 
ary phase. The impulsive (or transient) phase is a spiky, irregular 
signal, with multiple local extrema indicating: the arrival time of 
the shock-cusp at the null point (at t = 20.6 s), the formation of 
the first horizontal current sheet (at t = 35.7s, j z < 0), the forma- 
tion of the first vertical current sheet (at t = 59.1 s, j z > 0) and the 
formation of the second horizontal current sheet (at t = 90.0 s). 



The (current-sheet) signal is further contaminated due to addi- 
tion of small currents related to the propagation of shock waves 
across the null point. We define the end of the impulsive phase as 
the formation time of the second horizontal current sheet, which 
in our C = 1 simulation occurs at / = 90.0 s. 
After t = 90.0 s, the evolution of j z {0, 0, t) is much cleaner and 
closer to (damped) sinusoidal. For t > 90.0 s, the extrema exactly 
match the formation of the cyclic current sheets, with j z < 0/> 
indicating horizontal/vertical current sheets, respectively. We de- 
fine this 'cleaner signal' regime as the stationary phase, i.e. 
the regime characterised by a (damped) sinusoidal signal after 
the transients of the impulsive phase have leaked away, and this 
phase starts at the formation time of the second horizontal cur- 
rent sheet. The black dotted lines in Figure Q] indicate the forma- 
tion times of all the horizontal current sheets (for both phases). 
Note that in labelling these two regimes, i.e the impulsive phase 
and the stationary phase, we have adopted the terminology usu- 
ally associated with the excitation and damping of trapped and 
leaky modes in coronal loop oscillations (see, e.g., Terradas et 
al. 120051 120061 Luna et al. l2008t McLaughlin & Ofman|2008] 
and reference therein). 

3.1. Impulsive Phase 

Let us consider the impulsive phase of the evolution of j z (0, 0, f), 
i.e. evolution over < t < 90 (indicated by red line in Figure[T]l. 
We see that j ; (0, 0, t) is zero (since the equilibrium is potential) 
until the arrival of the fast oblique magnetic shock (which has its 
own associated current density) at the null point, indicated by the 
first (negative) extrema at t — 20.6s. This is followed by a second 
(local) minimum at t — 35.7 s, which is the formation time of the 
first horizontal current sheet (i.e. not at t — 20.6 s). The corre- 
spondence between extrema and current-sheet formation cannot 
be deduced from Figure Q] alone, but can be determined by com- 
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Fig. 2. Plot of (a) formation time (measured in seconds) and (b) length (measured in km and evaluated at the formation time) of 
first horizontal current sheet versus initial velocity amplitude (measured in km/s), where Vo = 25.2 km/s corresponds to C = 1 
simulation. 



parison with the evolution of the separatrices and contours of 
v. The value of j z (0,0,t = 35.7) = -41.25 j = -3.3 mA is 
proportional to the length of the first horizontal current sheet, 
which is measured directly from the numerical simulation as 
0.1 162 L = 1 16.2 km (for the C = 1 simulation). 
The formation time and length (at formation time) of the first 
horizontal current sheet is of key interest here, since this di- 
rectly reflects the driver (akin to a forcing term) of the sys- 
tem. We investigate how the formation time and length of the 
first horizontal current sheet vary as functions of the initial driv- 
ing amplitude (from equation[3]). Figure^ shows the formation 
time (measured in seconds) of the first horizontal current sheet 
as a function of initial velocity amplitude (measured in km/s), 
where Vo = 25.2 corresponds to C — 1 simulation. We see that 
the greater the driving amplitude, the earlier the first current 
sheet forms. This is intuitive as we expect fast oblique magnetic 
shocks with a larger amplitude to propagate faster and thus reach 
the null point more rapidly, compared to waves driven with a 
smaller amplitude. 

Figure |2]3 shows how the current sheet length (measured in km) 
evaluated at the formation time of the first horizontal current 
sheet varies as a function of initial velocity amplitude. We see 
that the greater the driving amplitude, the longer the horizontal 
current sheet. The length of the current sheet is also directly pro- 
portional to the value of j z (Q, 0), and so we also conclude that the 
greater the driving amplitude, the stronger the value of \j z {0, 0)| 
at the corresponding time. Again, this result is intuitive; it is the 
fast oblique magnetic shock that physically deforms the X-point 
into an osculating field structure (and ultimately into a horizon- 
tal current sheet) and thus we would expect stronger fast oblique 
magnetic shocks (i.e. with a larger amplitude) to deform, i.e. re- 
fract B away from the normal, and 'squash' the magnetic field 
to a greater extent, and thus to form longer and stronger current 
sheets. 

Finally, let us investigate the time taken to evolve from the first 
horizontal current sheet (at t = 35.7 s) to the formation time of 
the second horizontal current sheet (at t = 90.0 s for the C = 1 
simulation), namely the time taken for one complete cycle (i.e. 
horizontal current sheet evolves to vertical, evolves back to hor- 
izontal). This time, which we refer to as the impulsive period, is 
calculated to be t — 54.3 s (for initial amplitude vo = 25.2 km/s 
in the C — 1 system). Again, we now investigate how this im- 
pulsive period varies with the initial driving amplitude (oc Cvo) 
and this can be seen in Figure [3] Here, we see that the greater 
the initial driving amplitude, the shorter the resulting impulsive 
period, i.e. the shorter the time taken to evolve from the first 
horizontal current sheet to the second. Recalling the dependency 
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Fig. 3. Plot of impulsive period, i.e. time taken to evolve from 
the first horizontal current sheet to the second horizontal cur- 
rent sheet (measured in seconds) versus initial velocity ampli- 
tude (measured in km/s), where vo = 25.2 km/s corresponds to 
C — 1 simulation. 



seen Figure |2j), this means that for longer current sheets, the 
impulsive period is shorter. This means that the restoring force 
must be stronger and so we conclude that longer current sheets 
have a stronger restoring force. In this way, the system acts as a 
harmonic oscillator, i.e. the greater the displacement away from 
equilibrium, the stronger the restoring force. 

3.2. Stationary Phase 

Let us now investigate the stationary phase of the oscillation 
seen in Figure Q] (black line), i.e. the time evolution of j z (Q, 0, f) 
for t2 < t < 480 s, where ?2 is the formation time of the second 
horizontal current sheet {to = 90.0 s for the C — 1 simulation). 
FigureQ]also has an insert showing the (same) time evolution of 
j z (0, 0, t) for 158.9 < t < 480 s, i.e. same horizontal time axis, 
different vertical axis. Note that the insert does not show the start 
of the stationary phase, only a later part of it (starting at the time 
of the third horizontal current sheet). The green line indicates 
./final, i-e. the finite amount of current density left in the system 
at t - 480 s when the system has reached its final, non-potential 
state. For the C = 1 system, this final state occurs at t — 480 s (8 
mins) and j nna i = 0.8165 jo = 6.5 x 10~ 5 A. 
We see that there is clear oscillatory behaviour in the stationary 
phase and, moreover, the oscillation is exponentially decaying. 
This can be seen in Figure[T]and the blue dashed lines indicate an 
exponentially-damped envelope max( stationary x<?~' )f + ja n . d \ and 



5 



McLaughlin et al.: On the Periodicity of Oscillatory Reconnection 



(a) 
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initial amplitude [km/s] 



(b) 




0.04 

20 10 60 80 '00 120 
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Fig. 4. Plot of (a) stationary period, i.e. the time taken to evolve from the second horizontal current sheet to the third (measured in 
seconds) and (b) jsaai, i.e. j z (Q, 0) evaluated at t — 480 s, (measured in milliAmps) versus initial velocity amplitude (measured in 
km/s), where Vrj = 25.2 km/s corresponds to C = 1 simulation. 



mm Oz)lstationary x e At + y'g na i where A = -0.015 s 1 is determined 
experimentally. 

Let us now investigate the period associated with the station- 
ary phase, which we define as the time taken to evolve between 
horizontal current sheets. Specifically, we define the stationary 
period as the time taken to evolve from the second horizontal 
current sheet to the third, i.e. the first complete oscillation within 
the stationary phase. For the C — 1 simulation, these formation 
times are 90.0 s and 158.9 s resulting in a stationary period of 
69.0 s (note we present results here correct to 1 decimal places, 
but calculate periods to a greater degree of accuracy). Similar 
results are obtained for alternative definitions of the stationary 
periods, e.g. time taken to evolve from one vertical current sheet 
to the next. 

We now investigate how the stationary period varies with the 
initial driving amplitude and this can be seen in Figure^. Here, 
we see that the greater the initial driving amplitude, the shorter 
the resulting stationary period. Thus, these results are in agree- 
ment with those in {33.11 i.e. the stronger the initial driving am- 
plitude, the longer the resulting current sheet, thus the stronger 
the restoring force, thus the shorter the resulting period. Coupled 
with the exponential decay, we see that in the stationary phase, 
the system acts akin to a damped harmonic oscillator. 
We also measure all the proceeding oscillations in the stationary 
phase, i.e. time taken to evolve from the second/third horizontal 
current sheet to the third/fourth horizontal current sheet, and ob- 
tain similar periods of 67.8s and 66.6s respectively. Interestingly, 
this means that the stationary period appears to be slightly de- 
creasing by roughly 1.8% per oscillation. 
Finally, we investigate jfj n ai, tne finite amount of current density 
left in the system when the system has reached its final, non- 
potential state and this can be seen in Figure H]}. For the C — 1 
system, = 0.8165y'o = 6.5 x 10~ 5 A which is measured at 
t = 480 s. In all our simulations, > indicating that the 
(final) X-point is very slightly closed up in the vertical direction, 
i.e. j z > is associated with vertical current sheets. This is be- 
cause the (local) plasma to the left and right of the X-point is 
slightly hotter, since that is where the initial, strongest jet heat- 
ing occurs. Thus, the existence of this thermal-pressure gradient 
coupled with force balance requires the final state to be non- 
potential. 

From Figure |4j), we see that the stronger the initial driving am- 
plitude, the greater the value of jemi. Again, this is intuitive: fast 
oblique magnetic shocks with a larger amplitude will intersect 
to form stronger, hotter jets to the left and right of the X-point. 
Thus, this local plasma will be hotter at the end of the simula- 



tion, indicating a stronger thermal-pressure gradient and thus, in 
order to achieve force balance, a greater absolutely value of the 
Lorentz force, i.e. a greater value of ./'final- 



4. Conclusions 

This paper describes an investigation into the periodicity of os- 
cillatory reconnection, specifically oscillatory reconnection ini- 
tiated by a nonlinear fast magnetoacoustic wave deforming a 2D 
magnetic X-point. We have solved the compressible, resistive, 
nonlinear MHD equations using a Lagrangian remap, shock- 
capturing code (LARE2D) and have followed the numerical set- 
up of McLaughlin et al. (2009). As in that paper, we find that 
the fast magnetoacoustic wave develops into a fast oblique mag- 
netic shock wave which significantly deforms the local magnetic 
fieldlines, to the extent that the incoming wave deforms the null 
point into a cusp-like point which in turn collapses to a current 
sheet. The system then evolves periodically through a series of 
horizontal and vertical current sheets with associated changes in 
connectivity, i.e. the system demonstrates the mechanism of os- 
cillatory reconnection. 

The main focus of this paper is on the periodic nature of this 
oscillatory cycle of horizontal and vertical current sheets. For 
the first time, we identify two distinct phases in the oscillation: 
a transient, impulsive phase, encompassing the development of 
the first horizontal current sheet, the formation of the first ver- 
tical current sheet, and ending with the formation of the second 
horizontal current sheet. We define the stationary phase to begin 
at the formation of the second horizontal current sheet and thus 
this phase includes all the proceeding cyclic behaviour. 
In the impulsive phase, we find that greater the driving am- 
plitude (Cvo) of the velocity initial condition (equation [3), [a] 
the earlier the first horizontal current forms, [b] the longer its 
maximum length and [c] the greater its maximum current den- 
sity. These results are intuitive since we would expect magnetic 
shocks with larger amplitudes to propagate faster and thus arrive 
at the null point more rapidly than those with smaller amplitudes. 
We would also expect magnetic shocks with larger amplitudes to 
deform the pre-existing magnetic field to a greater extent (specif- 
ically to refract B away from the normal to a greater extent) and 
thus, ultimately, to form longer and stronger current sheets. 
We also investigate the time taken to evolve from the first hori- 
zontal current sheet to the second, which we labelled as the im- 
pulsive period. We find that the greater the initial velocity am- 
plitude (Cvo) the shorter the resultant impulsive period. Coupled 
with the results on current sheet length, this means that longer 



6 



McLaughlin et al.: On the Periodicity of Oscillatory Reconnection 



current sheets have shorter corresponding impulsive periods. In 
this way, the system is acts as a harmonic oscillator, i.e. the 
greater the displacement away from equilibrium, the stronger the 
restoring force, and thus the shorter the impulsive period. 
For a driving amplitude of 25.2 km/s (corresponding to C = 1 
simulation) we measure an impulsive period of 54.3 s. We also 
investigate the resultant impulsive periods for driving amplitudes 
6.3 - 126.2 km/s and find associated impulsive periods in the 
range 46.9 -65.4 s. 

The stationary phase is found to be dominated by an 
exponentially-decaying oscillation, tending to a finite value, 
./final- We also investigate the stationary period, namely the time 
taken to evolve from the second horizontal current sheet to the 
third, i.e. the first complete cycle of the stationary phase. As in 
the impulsive phase, we find that the greater the initial velocity 
amplitude (Cvo) the shorter the resultant stationary period. This 
is explained just as before, the greater the initial amplitude, the 
longer and stronger the current sheets at each stage, and thus the 
greater restoring force, leading to shorter periods (compared to 
smaller initial amplitude, shorter resultant current sheets, weaker 
restoring force and thus longer periods). Hence, again the system 
acts as a harmonic oscillator and, coupled with the exponential 
decay, we relate the oscillatory reconnection mechanism to that 
of a damped harmonic oscillator during the stationary phase. 
For a driving amplitude of 25.2 km/s (corresponding to C — 
1 simulation) we measure a stationary period of 69.0 s. We 
also investigate the resultant stationary periods for driving am- 
plitudes 6.3 - 126.2 km/s and find associated stationary peri- 
ods in the range 56.3 - 78.9 s, i.e. these are high frequency 
(0.0127 - 0.0178 Hz) oscillations. 

It is also prudent at this stage to ask what determines this sta- 
tionary period and what determines the exponentially-decaying 
timescale. To this end, we can consider the work of Craig & 
McClymont (119911 1 who investigated the relaxation of a 2D X- 
point disturbed from equilibrium. By neglecting both nonlinear 
and thermal pressure effects, Craig & McClymont ( [19911 1 de- 
rived an analytical prediction for two timescales: 

^oscillation ~ 2 In R m , fdecay ~ ^oscillation 

where we identify filiation as our stationary period and fdecay as 
our decay time; I /A. For a driving amplitude of 25.2 km/s, these 
correspond to filiation ~ 109.6 s, compared to our measured 
stationary period of 69.0 s, and fdecay ~ 76.7 s compared to our 
measured decay time of 1/A = 1/0.015 = 66.7 s, given that in 
our investigations R m = 10 4 and time is made dimensional using 
to — 7.93 s. Thus, given the (relative) simplicity of the Craig & 
McClymont (119911 1 system, these estimates are in fair agreement 
with our results. Note, however, that these simple analytical for- 
mulae cannot predict the variation in period versus amplitude 
of the initial velocity driver. This suggests that nonlinear effects 
and thermal-pressure gradients play a crucial role, which seems 
reasonable given that the restoring force of oscillatory reconnec- 
tion has been shown to be a dynamic competition between the 
thermal-pressure gradients and the Lorentz force (see §3.2 of 
Murray et al. 120091 §3.3 of McLaughlin et al. 120091 Figure 7 of 
Threlfalletal.|2pT2j. 

It is also important to note that there is a significant difference 
between the impulsive period (e.g. 54.3 s for C — 1 system) and 
the stationary period (e.g. 69.0 s for the C = 1 system), and that 
in every numerical experiment we find that the stationary period 
is longer than the impulsive period. This indicates that different 
physical processes dominate in each phase (i.e. the deformation 
of the X-point by the shock and the importance of jet heating 
in the impulsive phase, and the elastic motion of the magnetic 



field trying to get back to equilibrium in the stationary phase) 
and validates our approach of dividing the whole evolution into 
two distinct phases. 

This difference in impulsive and stationary periodicities actually 
has an intriguing caveat : it is important to note that when oscil- 
latory reconnection is seen in, say, a numerical simulation, one 
must be careful to interpret which phase one is actually observ- 
ing and to observe several oscillations, i.e. if only two periods 
are seen, say the impulsive period followed by a single station- 
ary period, then one would conclude that the period was actually 
increasing between oscillations. A similar result would pertain 
in solar observations of oscillatory reconnection, i.e. the first pe- 
riod measured would be shorter than the proceding periods (as- 
suming the first, i.e. impulsive, period is also observed). This is 
a clear prediction for the oscillatory reconnection mechanism. 
It is also important to note that, as shown by McLaughlin et al. 
(120121 1. the mechanism periodically generates v t and v v but that 
these are generated exponentially damped. Thus, if such signals 
are detected, then they may be decaying not due to a particu- 
lar damping mechanism, but due to the generation mechanism 
itself. 

In addition to the stationary period, we also measured all the pro- 
ceeding periods in the stationary phase, e.g. time taken to evolve 
from third horizontal current sheet to fourth, etc. Interestingly, 
it was found that the period very slightly decreases by roughly 
1.8% per oscillation. The exact reason for this decrease in the 
period is uncertain and may be a numerical effect. This will be 
investigated further in future work. 

As in McLaughlin et al. (2009), it was found that the final state 
(i.e. velocity zero, oscillatory behaviour ceased) is in force bal- 
ance but is non-potential and a small, finite amount of current 
density exists in the system, /final. The (final) X-point is very 
slightly closed up in the vertical direction, i.e. j z > is as- 
sociated with vertical current sheets. This is because the (lo- 
cal) plasma to the left and right of the X-point is slightly hot- 
ter, since that is where the initial, strongest, jet heating occurs. 
Thus, the existence of this thermal-pressure gradient in force bal- 
ance requires the final state to be non-potential. We find that the 
greater the initial velocity amplitude (Cvo) the larger the value 
of /final- Again, this is intuitive: fast oblique magnetic shocks 
with a greater amplitude will overlap to form stronger, hotter 
jets to the left and right of the (equilibrium) X-point. Thus, this 
local plasma will be hotter at the end of the simulation, indi- 
cating a stronger thermal-pressure gradient and thus, in order to 
achieve force balance, a greater value of the Lorentz force, i.e. a 
larger value of j 6ma . 

We have presented an investigation into the periodic nature of 
oscillatory reconnection and have found that an aperiodic driver 
can naturally generate a period signal via the physical mecha- 
nism of oscillatory reconnection. We have found that the system 
behaves akin to a damped harmonic oscillator. Again, this is 
not surprising: effectively our velocity initial condition can be 
thought of as injecting a finite amount of energy into the os- 
cillatory reconnection mechanism and so intuitively the resul- 
tant periodic behaviour must be finite in duration, i.e. this is a 
dynamic reconnection phenomena as opposed to the classical 
steady-state, time-independent reconnection models. 
Oscillatory Reconnection may also play a role in generating 
quasi-periodic pulsations (see, e.g., reviews by Aschwanden 
2003; Nakariakov & Melnikov 2009). Oscillatory behavior has 
been reported in a number of solar and stellar flare observations 
(e.g. Mathioudakis et al. [20031 [20061 M cAteer et al. l2005l Inglis 
et al. 120081 Inglis & Na kariako v 120091 Nakaria kov et al. l20T01 
Nakariakov & Zimovets I2PTT1 Inglis & Dennis 12PT21 Shen & 
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Liu l2012t but the generation mechanism responsible remains an 
open question. 

We believe the physical mechanism of oscillatory reconnec- 
tion described in this paper is a robust, general phenomenon 
that will be observed in other systems that demonstrate finite- 
duration/non-steady-state reconnection (although we have only 
presented a specific example of oscillatory reconnection in this 
paper). For example, evidence of oscillatory reconnection in 3D 
flux emergence simulations has been reported by Archontis et al. 
(I20T01 I. 
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